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Abstract 

In this paper we propose to interpret the large discretization arti- 
facts affecting the neutral pion mass in maximally twisted lattice QCD 
simulations as 0(a 2 ) effects whose magnitude is roughly proportional 
to the modulus square of the (continuum) matrix element of the pseu- 
doscalar density operator between vacuum and one-pion state. The 
numerical size of this quantity is determined by the dynamical mech- 
anism of spontaneous chiral symmetry breaking and turns out to be 
substantially larger than its natural magnitude set by the value of 
Aqcd- 



1 Introduction 



In recent simulations of maximally twisted lattice QCD (Mtm-LQCD) pQ 
with Nf = 2 dynamical light quarks, numerical results [2, [3] on the neutral 
pseudoscalar meson mass at different lattice spacings (in the range from 0.10 
to 0.065 fm) and quark masses (corresponding to charged pseudoscalar meson 
masses ranging from 600 to 300 MeV) show large cutoff effects. The latter 
appear at fixed lattice spacing in the form of sizeable deviations from the 
expected continuum QCD chiral behaviour. 

This finding is in striking contrast with the smallness of cutoff effects 
observed not only in the mass of the charged pions (which is related through 
the Goldstone theorem to exactly conserved lattice currents [1]), but also in 
all the other hadronic observables so far measured by the European Twisted 
Mass (ETM) Collaboration (see sect. [3] and refs. El El El El GO] ) ■ Quite 
remarkably small lattice artifacts are found even in matrix elements where 
the neutral pion is involved (see Table [T] below) . 

In Mtm-LQCD the additive cutoff effect in the squared mass of the neu- 
tral pion is expected on general grounds [H] to be of the order a 2 AQ CD , 
whereas the leading additive corrections to the squared mass of the charged 
pion are [121 [El H] only of order a 4 Aq CD and a 2 /i g AQ CD (fi q denotes the 
bare quark mass). The neutral vs. charged pion (squared) mass splitting, 
Am 2 | ^ tm = 777.^3 1 1 — m^ ± | l, in the small quark mass region is then dom- 
inated by 0(a 2 AQ CD ) terms. Numerical data (see Fig. [T] and Table H]) at 
two different lattice spacings and for charged pion masses below 500 MeV 
are indeed consistent with the expected behaviour of AmJ|f tm , but with a 
large coefficient in front of the 0(a 2 AQ CD ) correction 0. For instance, taking 
Aqcd = 250 MeV, one finds Am 2 |^ tm 50<2 2 Aq CD . 

On a more general ground, the question of the numerical size of the 
neutral vs. charged pion mass splitting is an important issue to assess the 
viability of the Mtm-LQCD approach. This is so mainly because of the 
established relation [TJl [TJ] between the magnitude of the pion mass splitting 
and the strength of the metastabilities detected in the theory at coarse lattice 
spacings EH El LIE]. 

One may suspect that the size of Am 2 |^ tm represents the generic magni- 
tude of the isospin breaking effects inherent in the twisted form of the action. 
This is not so, however. Numerical data in several other observables indi- 
cate, in fact, that in general isospin breaking artifacts are pretty small, with 
a relative magnitude in line with their naive order of magnitude estimate. 

1 The issue of the Nf dependence of this coefficient is beyond the scope of this paper 
and will be discussed elsewhere. 
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The main purpose of this paper is to provide an explanation for this 
peculiar pattern of isospin breaking effects, and in particular for the fact 
that only the neutral vs. charged lattice pion mass splitting among all the 
other 0(a 2 ) effects is large. 

Relying on arguments based on the Symanzik analysis [US, [20] of lattice 
artifacts, we are able to identify the form of the leading 0(a 2 ) corrections 
that in Mtm-LQCD affect the value of the squared mass of the neutral pseu- 
doscalar meson. Interestingly one can give an approximate evaluation of 
these effects finding that they are proportional to the modulus square of the 
continuum matrix element of the isotriplet pseudoscalar density operator be- 
tween vacuum and one-pion state. The latter is a dynamically large number 
determined by the mechanism of spontaneous chiral symmetry breaking. 

Lattice artifacts of course also depend on the many unknown coefficients 
multiplying the operators occurring in the Symanzik expansion. However, 
the numerical evidence provided in this paper (see sect. [3]) and in refs. [3J, [6] 
about the fact that all the other so far measured observables exhibit only 
small lattice artifacts should be taken as a strong indication that the coeffi- 
cients multiplying the (quark-mass independent) operators of the Symanzik 
effective Lagrangian of relevance here (operators of dimension 5 and 6) are 
not unnaturally large, at least in the gauge coupling regime explored up to 
now in ETMC simulations. Indeed, were not this the case, one would have 
seen large cutoff effects in some other physical observables besides the neutral 
pion mass. 

A preliminary version of this investigation was presented some time ago 
in [2T]. Ideas and a few results along the line of reasoning developed in this 
paper have been already put forward in j22j [121 02] and in [23] . 

1.1 Plan of the paper 

The plan of the paper is as follows. In sect. [2] we illustrate the nature and 
the structure of 0(a 2 ) artifacts in Mtm-LQCD and how they affect charged 
and neutral pion masses. In sect. [3] we give numerical evidence that among 
the many observables recently measured by the ETM Collaboration only 
the pion mass splitting seems to display large cutoff effects. Conclusions 
can be found in sect. [H More technical issues are discussed in Appendices. 
In Appendix A we review the notion and the properties of the Symanzik 
expansion for the description of the cutoff effects of LQCD with Wilson 
fermions (either maximally twisted or untwisted) and how, depending on the 
value of the twisting angle, its form is affected by the way the critical mass is 
determined. In Appendix B we give the structure of the Symanzik expansion 
of the two-point pseudoscalar correlator from which the formula for the 0(a 2 ) 
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discretization errors affecting the pion masses can be derived. Finally in the 
(long) Appendix C we give details on the theoretical analysis and numerical 
estimate of the matrix elements controlling the magnitude of certain cutoff 
corrections affecting neutral and charged pion masses. 



2 Neutral and charged pion mass in Mtm-LQCD 

In this section we want to describe and compare the structure of the lattice 
artifacts affecting the neutral and the charged pion mass in Mtm-LQCD. 
The key formulae yielding the 0(a 2 ) corrections to the charged and neutral 
squared pion mass are immediately derived by writing down the leading cor- 
rection to the pion (rest) energy induced by the 0(a 2 ) terms of the Symanzik 
expansion. In Mtm-LQCD one gets (b = ±, 3) 



2 , 2 

m T + a 



(vr b (0)|£f m (0)|7r b (0)) 



+ 



;(7T b (0)| / d 4 xCf l m (x)£f m (0)|7T 6 (0')) 



+ 0(a 2 m 2 a 4 ) 



(2.1) 



where the form of £j^ tm and £^ tm is detailed in Appendix A. Eq. ( 12. ip com- 
plies with elementary lowest order perturbation theory formulae and Lorentz 
covariance. For completeness in Appendix B we give another derivation of 
eq. (12. ip starting from the Symanzik expansion of the Fourier Transform of 
the two-point pseudoscalar correlator (no sum over 6), 



ri(p) 



a ^e ipx {P b (x)P b (0)) 



(2.2) 



where P b = ^z7s(r fc /2)^ and t\ b = 1,2,3, are the Pauli matrices. 

The main conclusions of this paper are derived from the following two 
observations concerning the structure of eqs. (12.11) . 

1) The pion matrix elements (vr b (0)|£^ tm (0)|7r 6 (0)) | cont are 0(1) for 6 = 3, 
but only 0(m 2 ) for b — ± (see Appendix B), viz. 



(7^(0)14 



7^(0)) 



0(m 



C. = (7r 3 (o)|/:f m (o)k a (o)) 



cont 



nJ 1 

0(1) 



(2.3) 
(2.4) 



2) Upon use of the optimal critical mass [3], one can show (see Ap- 
pendix C) that the matrix element 



^55 



d 4 xCf tm (x)Cf tm 



Ao)) 



(2.5) 
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is parametrically of 0(m 2 ) in the chiral limit for b — ±, and numerically 
negligible (in modulus) compared to \( n \ in the case of b = 3. 

The direct consequence of these statements is that close to the chiral limit 
there are no additive 0(a 2 ) terms affecting the value of the lattice charged 
squared pion mass. The neutral squared pion mass instead receives non- 
vanishing corrections at this order. In formulae one gets 



L = ml + a 2 A% + 0(a 2 ml, a 4 ) = ml + 0(a 2 m 2 , a 4 ) , (2.6) 

,2 i „2l> I a 3 i i r\( Jl^l 



l = m 2 +a% + A^ + 0(a 2 m 2 ,a 4 )- (2.7) 

One can also show that in Ag 5 the numerically dominant 0(m 2 ) terms, com- 
ing from the insertion of one-pion states in (12. 5p . are actually independent 
of the isospin index b. As a result what is left in the difference Ajj 5 — Af 5 is 
a tiny correction, so that for the pion squared mass difference we can finally 
write down the simple formula 



Am 2 



Mtm „ 






= "0 


— mti 




L 77 


L n 


L ~~ 



a 2 ( n + 0(a 2 ml,a 4 ) , (2. 



where the symbol ~ is to remind that the (negligibly small) A^ 5 — A^ 5 cor- 
rection has been dropped. 

In the rest of the paper we want to give a proof of the propositions 1) and 
2) above and provide a numerical estimate of A55 and ( n . In this section we 
shall see in particular that for dynamical reasons ( n is large compared to its 
natural value (~ Aq CD ), providing in this way an explanation for the fairly 
big value of Am 2 \^ tm measured in numerical simulations [2, 3J. 

The arguments leading to the eqs. (I2.6p - fl2.7l) and the statements above 
are quite elaborated. To keep the line of reasoning as straight as possible, we 
have deferred them to Appendix B (derivation of eqs. (12.61) and (12.71) ) and 
Appendix C (expression and estimate of A^). 



2.1 Estimating 0(a 2 ) lattice artifacts in Am 2 ^ 

Relying on the results of Appendix C about the numerical irrelevance of A^ 5 
(either because parametrically of 0(m 2 ), for b = ±, or because numerically 
small, for b = 3), we only need to evaluate ( n in (the chiral limit of) con- 
tinuum QCD in order to estimate the size of the 0(a 2 ) artifacts in ( 12. 7p 
and (12. 8p . This can be done under the assumption that a sufficiently accu- 
rate order of magnitude estimate of this hadronic parameter can be obtained 
in the vacuum saturation approximation (VSA). Quenched LQCD studies 
give support to the assumption that VSA works well for matrix elements of 
four-fermion operators between pseudoscalar states [2U[25]. It is very likely 
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that this remains true also in the unquenched theory. Indications in this 
sense are actually born out by recent studies with two or three dynamical 
flavours pEJETj. 

2.1.1 The theoretical argument 

The evaluation of ( n (see eq. (12.41) ) requires the estimate of the matrix ele- 
ments of the operators in eq. (1A.8I) (rotated from the twisted to the physical 
quark basis) between zero three-momentum neutral pion states. This is done 
in two steps. 

I) By the use of classical soft pion theorem's (SPT's) [2H1 EH] one rec- 
ognizes that some operators have 0(1) neutral pion matrix elements, while 
others vanish proportionally to m^. The latter are those that are invariant 
under axial-r 3 transformations. The former are 



The operators (I2.9P are obtained from those of eq. (IA.9I) after use of the chi- 
ral rotation ( 1A.2I) . We note that SPT's relate through the trivial numerical 
factor —1 the matrix elements between neutral single-pion states of the two 
operators in each line H. Not surprisingly this remark implies that would 
vanish if we were to deal with an exactly chiral invariant lattice formulation. 
In fact, such a formulation would be described by a Symanzik effective La- 
grangian where the two operators in each line of the list (12. 9p . being related by 
a chiral transformation, would necessarily appear with identical coefficients. 

II) In VSA one can give an estimate of the value of (eq. (12.41) ) as 
follows. First of all, one has to rewrite the matrix elements of the opera- 
tors (12.91) between zero-momentum neutral pion states in the form that is 
obtained reducing external pions by the use of SPT's H In this way one 

2 For the reader convenience we recall that SPT's amount to the relation i/,r(a + 
n b (0)\O c \(3) = (a\[Q b A , O c ]\(3), where Q b A is the axial charge with isospin index b. External 
states must carry isospin indices such that the two matrix elements are not trivially van- 
ishing. We also note the crossing symmetry relation (a + n b (0)\O c \j3) = (a\O c \/3 + TT b (0)). 
For an introduction to SPT's and a more detailed discussion (including limitations in 
their use if intermediate states with the quantum numbers of O c \a), or O c \(3) happen to 
be degenerate in energy with the states |a) or |/3)) see e.g. [30], sect. IV-5, and references 
therein. 

In doing so we are enforcing consistency with SPT's in using VSA-based estimates 
of the matrix elements of the operators (|2.9p . This step is necessary because it is not 
guaranteed that VSA estimates comply with SPT relations. 



I (^75^0 (V'*75J0) , 



(2.9) 
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checks that only the pion matrix elements of the two operators in the first 
line of the list (12.91) are not vanishing in VSA. Their actual contribution to 
C,r depends of course on the magnitude of the Symanzik coefficients (dimen- 
sionless functions of the bare gauge coupling) by which the two operators are 
multiplied in the expression of £j^ tm . In the gauge coupling regime of interest 
for large volume simulations of Mtm-LQCD, these coefficients are expected 
to be 0(1) quantities, but their value (and sign) is otherwise unknown. We 
must thus limit ourselves to an order of magnitude estimate of for which, 
relying on the line of reasoning outlined above, we write 

ICI^ A ^-2^, (2.10) 

where ~ 92.4 MeV. The quantities E and G n denote the continuum renor- 
malized chiral condensate and the vacuum-to-pion matrix element of the 
pseudoscalar density, with 

P 6 = V^7 5 tV>, (2-11) 
S° = ^. (2.12) 

The first relation, as explicitly indicated, comes from a direct use of SPT's 
and vacuum insertion. The second holds up to 0(m\) corrections and comes 
from combining the Gell-Mann-Oakes-Renner formula [31] 

2/2 g E = flml + 0(m*) (2.13) 

with the continuum Ward-Takahashi identity (WTI) 

2fi q G^ = Uml. (2.14) 

We conclude this subsection by stressing that sign of ( w , which, as we said, 
depends on the operator coefficients in £^ tm , is not predicted by our argu- 
ments and is to be learned from numerical simulation experience. 



G n = (^|PV(0)> 

Cl 

X 1 2 cont ' 



2.1.2 Numerics for eq. ^2. 1 (A) 

A numerical evaluation of G^, or equivalently of can be obtained in 

various ways, i.e. either exploiting the direct lattice measurement of G n or 
by using (at the physical pion point) the continuum WTI f)2.14p . 

I) The direct lattice measurements of carried out in ref. [21 [3], [9] at 
different lattice resolutions give in the continuum and chiral limit 

GV(MS, 2 GeV) ~ (490 MeV) 2 , [-t(MS, 2 GeV)] 1/3 ~ 275 MeV (2.15) 
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with a total error on G n not larger than 5%. In order to arrive at these 
numbers the renormalization constant of the operator P b , computed in the 
RI-MOM scheme [32] and converted to the MS scheme, was employed [33J. 

II) A second, independent determination of G n can be obtained by making 
reference to only continuum quantities if eq. (12.141) and the PDG [M] estimate 
jl q (MS, 2 GeV) = [{rh u + m d )/2}(MS, 2 GeV) ~ 3.8 MeV are used. In this 
way one gets at the physical pion mass 



G„(MS, 2 GeV) ~ (470 MeV) . (2.16) 

Taking into account the small error in the extrapolation from the chiral point 
to the physical pion and the uncertainty on the light quark mass value, we 
see that the two estimates of G n are in very good agreement with each other. 

The conclusion of this analysis is that numerically the estimate (12.81) of 
the neutral to charged squared mass shift is substantially larger than its 
natural size, a 2 AQ CD , by a factor (see eqs. (I2.10p . (12.151) and (12.161) ) of the 
order of 2G 2 /A^ CD ~ 25 30. 

2.2 Numerical results for Am 2 |r 

In this section we wish to summarize the simulation results for Am 2 ^ that 
have been obtained from Nf = 2 ETMC ensembles at two different lattice 
resolutions, a -1 ~ 2.3 GeV and a -1 ~ 2.8 GeV (corresponding to j3 = 3.9 and 
(3 = 4.05, respectively (see sect. IG.3.21 and Table [21 for details), and compare 
them with the number expected on the basis of eq. (I2.10p for ( v as well as 
with the estimate of Ag 5 — A 55 derived in Appendix C (see eq. (1C.26I) ). 

In Fig. CE] we display the lattice results for (r /a) 2 rQAm 2 \l vs. TqIji^Il 
for the six gauge configuration ensembles specified in Table [2] The factor of 
(r /a) 2 in the vertical axis has been inserted @ to be able to combine together 
data from f3 = 3.9 and j3 = 4.05, thus giving an impression of the quality of 
the scaling behaviour of Am 2 | l ■ In the same figure we also report the value of 
VqCtt (shaded square), which from eqs. (12.101) . (12.151) . (12.16!) an d r o — 0.44 fm 
is computed to be[f| 

4( n ~ -(2.4-2.8) @//, = 0. (2.17) 

The uncertainty in (I2.17p . which is dominated by the approximations in 
the theoretical arguments given in sect. 12.1.11 corresponds to the size of the 



According to ref. 0, we use r /a = 5.22(2) at (3 = 3.9 and r a /a = 6.61(3) at j3 = 4.05. 
5 The sign of (jr, which appears in eq. (|2.17p has been inferred from the numerical 
evidence that m^3 < m T ± . 
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Figure 1: Data for [r^j a^r^Am 2 ^ vs. ?"o^±|l are compared with the esti- 
mate of /"oCtf (shaded square) and ro(A| 5 — A^ 5 ) (green dot). 



symbol used in the figure. Finally the green dot in Fig. [T] represents the result 
of the estimate (IC.26I) we give of Af 5 — Af 5 expressed in r -units, yielding 

r\(A\ h - Af 5 ) = -0.10 ± 0.06 ± 0.03 , (2.18) 

where the first error is essentially statistical and stems from the evaluation 
of £ n in eq. (IC54[) and the second is the systematic uncertainty reflecting the 



approximations involved in the way we parametrize (sect. IC.ll) and evaluate 
(sect. EL2D the quantity A| 5 - A^. 

We see from Fig. [1] that within statistical errors (mainly coming from 
the noisy quark-disconnected diagrams entering the calculation of 7?v0 the 
available data for Am 2 \ l show the expected scaling behaviour with the lat- 
tice spacing. The quark mass dependence of (r /a) 2 rpAm^ \l turns out to 
be of reasonable magnitude and not inconsistent with our theoretical esti- 
mate (12.17ft within the large uncertainty affecting (dAm%/d£L g )\L (from the 
plotted data one finds (dAm 2 n /djl q )\ L ~ 300 900 MeV at /3 = 3.9). 
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The estimated value of Tq(A| 5 — A^ 5 ) with its error (see eq. (I2.18P ) is 
(in modulus) much smaller than the observed pion (squared) mass splitting, 
and also about twenty times smaller than ^"o I Ctt I • We thus conclude that 
the dominant contribution to comes from the a 2 ( n correction, as 

announced. 

3 0(a 2 ) isospin violating effects in Mtm-LQCD 

In this section we want to provide numerical evidence of the fact that in 
Mtm-LQCD simulations 0(a 2 ) isospin breaking effects are negligible in all 
the physical quantities measured up to now by the ETM Collaboration, with 
the exception of the neutral pion mass. 

In Tabled] we report the relative difference of pseudoscalar meson squared 
masses and decay constants between neutral and charged particles for pseu- 
doscalar and vector mesons, as well as the relative splitting between the 
masses of A + and A ++ baryons. Most of these data already appeared in the 
conference contributions of ref. {3j and in ref. [6]. 

We refer to [7] for a detailed description of the techniques used to com- 
pute correlators, extract physical observables and perform the necessary error 
analysis. The results presented in this paper are based on the configurations 
generated by the ETM Collaboration using the formulation of lattice QCD at 
maximal twist (with csw = |35j) and tree- level Symanzik improved gluon 
action. We show data coming from (3 = 3.9 [2] and (3 = 4.05 0, E] simula- 
tions at two, roughly matched (in physical units) values of the twisted mass 
parameter fi q for each (3. Again we stress that the computation of the mass 
and decay constant of neutral (pseudoscalar and vector) mesons involves the 
evaluation of quark-disconnected diagrams, which is the reason for the rela- 
tively larger errors affecting these quantities. In the case of A-baryon masses 
(where no computation of quark-disconnected diagrams is required) we give 
results at four values of fi g for (3 = 3.9 and (3 = 4.05. 

We quote in Table [T] the relative isospin splitting, R s , measured between 
the observables O and O', namely 

O-O' 

with O (O') denoting the quantity in the charged (neutral) sector in the case 
of meson data and the A +,0 -mass (A ++ '~-mass) in the case of nucleons. We 
recall that, owing to the invariance of the Mtm-LQCD action under parity 
and u <-> d flavour exchange, A + and A baryons, as well as A ++ and the 
A~, are exactly mass-degenerate. 
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It should also be remarked that to calculate i£ s [0] in the cases of the 
pseudoscalar and vector (isotriplet) meson decay constants use has been made 
of the appropriate renormalization constants (Za and Zy®) which were taken 
from ref. [33]. The expressions of the pseudoscalar and vector meson decay 
constants in the physical quark basis read (no sum over b) 



f psb = -±-(Q\^ l0l5 lr b ^) R \PS b ) , b = ±,3, (3.4) 
/ v * = -^-mi, lk \r b ^) R \V b ) , b = ±, 3 , (3.5) 



where the suffix R denotes renormalized operators and k is a spatial index Q. 
Preliminary results on fyb have appeared in [27] . 

The striking conclusion which emerges looking at Table [T] is that among 
the many relative splittings reported there we observe a large and statistically 
significant value of R s only when pion masses are compared. For all other 
observables the splitting is consistent with zero within the quoted errors. 



4 Concluding remarks 

In this paper we have proposed a theoretical explanation for the size of 
the 0(a 2 ) flavour violating cutoff effects responsible for the large splitting 
between neutral (m^a) and charged (m^±) pion mass seen in recent Mtm- 
LQCD dynamical simulations [2J. We have identified the origin of the ad- 
ditive 0(a 2 ) cutoff effects appearing in Trit 3 , which are absent in m^± as 
the latter has discretization errors only of order a 2 n q , a 4 and higher. The 
magnitude of the 0(a 2 ) lattice artifacts in the neutral pion mass is con- 
trolled, in the language of the Symanzik expansion, by the continuum matrix 
element (n 3 (6)\[Cf tm (0) - (1/2) J d A xCf tm (x)Cf tm (0)]\n 3 (6))\ cont , which, if 
Mtm-LQCD is defined via the optimal critical mass estimate, is numerically 

6 We recall that Za and Zy are the renormalization constants of the x-basis bare 
operators X7m75 t& X an d Xl^^Xi respectively. For the reader convenience we also report 
the renormalization rules for the currents in Mtm-LQCD [TTJ, in the physical basis. 
For the charged (b = ±) axial and vector currents they read 

(^7p75T ± V')i?. = ZV^7p75T ±- , (^7 AI r ± f/')fl = Za^^iI) , (3.2) 

while for their neutral (b — 3) counterparts one has 

(tpJt^l5T 3 tp)R = Z A 4>-y^5T 3 ip , (i/HVT 3 -!/;)/?, = Zyipj^T 3 ^ . (3.3) 



7 To be consistent with the H(4) invariance of the meson correlators the matrix elements 
involving charged (b — ±) or neutral (b = 3) mesons are normalized in eqs. (I3.4[) and (|3 . 5[) 
with the corresponding charged or neutral meson mass. 
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dominated by ( w = (7r 3 (0)|£^ tm |7r 3 (0))| cont . In the vacuum saturation ap- 
proximation one gets the theoretical estimate \( n \ ~ 1G\ ~ (560-^580 MeV) 4 , 
up to an unknown multiplicative factor depending on the details of the lattice 
action and the number of dynamical quark flavours (actually one can argue 
that this coefficient grows linearly with Nf). 

We have shown that numerical simulations of Mtm-LQCD with Nf = 2 
light dynamical quarks yield values of (r / ' a 2 ) 2 r^ j /S.m 2 K \L with negative sign 
and a magnitude that is in line with the above theoretical estimate. We have 
also provided evidence that among the many physical quantities recently 
measured in numerical simulations only the neutral pion mass appears to be 
affected by large flavour-breaking discretization errors. 

Given the independence on the twist angle of the Symanzik effective ac- 
tion in the limit of vanishing quark mass, these results may have a bear- 
ing also on LQCD simulations employing standard (i.e. untwisted) Wilson 
fermions. Indeed also in this formulation 0(a 2 ) cutoff effects of the same 
nature as those observed in Mtm-LQCD are expected to be present. In 
particular a formula completely analogous to eq. (12.11) holds true for the 
squared pion mass. Typically these 0(a 2 ) cutoff effects are reabsorbed in the 
definition of the current quark mass (which by construction vanishes simul- 
taneously with the pion mass) and thus in the corresponding definition of the 
critical mass. The consequences of this way of proceeding are in principle 
detectable in the form of 0(a 2 ) cutoff effects in the masses of all the other 
hadrons and, more generally, in all the quantities that have a non-negligibly 
small dependence on the quark mass. Actually, in what observables exactly 
they will appear depends on the precise way the critical mass happens to 
have been determined (see the analysis in Appendix A). 



Appendix A - Symanzik expansion and critical mass in 
LQCD with Wilson fermions 

We consider Nf = 2 LQCD with quarks regularized as Wilson fermions. For 
generic values of the bare (twisted, /j, q , and untwisted, m ) mass parameters 
the lattice action reads 



S, 



X 



~ a la q 

7 • V - - V V + csw—a ■ F + m + i/i ff 7 5 r- 



X 



(A.l) 



where V M = |[V M + V*], with V M (V*) the forward (backward) gauge co- 
variant lattice derivative. For the sake of generality we have also allowed for 



8 We adopt the conventions and notations of refs. [TTJ [33 [37J . In particular the quark 
fields in the "physical" and "twisted" basis will be denoted by ip, 4> and \, x, respectively. 
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the presence of the clover term. We are especially interested in two specific 
regularizations of the Dirac- Wilson action comprised in (1A.1D . 

The first is the Mtm-LQCD regularization which is obtained from I 1A.1I) 



by setting fi q = 0(a°) and m = M c e r , where M c e r is some estimate of the 
critical mass. The physical interpretation of this scheme is most transparent 
in the quark basis resulting from the field transformation 

if; = exp(i7T7 5 r 3 /4)x, $ = Xexp(i7T7 5 r 3 /4) , (A.2) 
where the lattice action takes the form 



7" V - i^T* f-^V*V + C SW JCT-F + M e cl ) + fi q 



^.(A.3) 



This quark basis is referred to as the "physical" basis, because the quark 
mass term takes the form n q ipip with \x q real [TTj . 

The second is the (clover) standard Wilson fermion action, S£', which is 
obtained by setting fi q = and mo = m + M£ with m an 0(a°) quantity. 
With this choice the most appropriate basis for discussing physics is the 
X-basis itself in which eq. ( lA.ip was written in the first place. 

A.l Symanzik expansion of Wilson fermion LQCD 



The Symanzik effective Lagrangian associated to the Wilson action I 1A.1I) 
reads 

A = C YM + X [P+ m + H 5 T 3 fx q }x , (A.5) 

5C Sym = aC 5 + a 2 C 6 + 0(a 3 ) , (A.6) 

where the four-dimensional operator, £ 4 , specifies the formal target theory 
in which continuum correlators are evaluated. The very definition of effec- 
tive action (all the necessary logarithmic factors are understood) as a tool 
to describe the a dependence of lattice correlators implies that the mass pa- 
rameters in £4, if not exactly vanishing, must be 0(a°) quantities. Thus 
all the lattice artifacts affecting M£ T will be described by higher dimensional 
operators in ££s y m of the form a fc £ fc AqcDXX> k > 1, with 5k dimensionless 
coefficients, which are functions of the gauge coupling. 

After using the equations of motion entailed by £4, the 0(a) piece of 
oXsym reads [201 13] 



C 5 = b, 



<5;SwXiv ■ F x + 5iAq CD xx + 0(m, fi q ) . (A.7) 
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The terms multiplied by powers of m and/or fi q are not specified in eq. (IA.7[) 



as they are not of relevance for the topic discussed in this paper. We recall 
that the coefficients &5 ; sw and 8± vanish if csw in eq. ( lA.lj) is set to the value 
appropriate for Symanzik 0(a) improvement. 

The 0(a 2 ) part of oXs y m has a more complicated expression, of the type 



id 



£ 6 = ]T Ki^iT + K^Xl^D^fx + Y, &6;i$6;i + £ 2 AQCDXX + 0(m, /i ff ) , (A.8) 



i=l i=5 



$6;5 = 


f(xxXxx) , 


$6;6 = 


= zE b (xT b x)(xr b x) , 

= iEfe(XH5^x)(X75^ 6 x) , 
= iEfc(X7AT b x)(X7AT fe x) , 

= lEb(XlXl5T b x)(XlXl5T b x) 
= \ Y,b(XVXvT b x)(X°XvT b x) ■ 


$6;7 = 




$6;8 = 


$6;9 = 


j(xi\x)(xi\x) , 


$6;10 


$6;11 = 


= i(^7a75X)(X7a7sX) , 


$6;12 


$6;13 = 


= \(WxuX)(W\uX) , 


$6;14 



where the first three operators are purely gluonic, the fourth is a non-Lorentz 
invariant fermionic bilinear and the remaining ones are four-fermion opera- 
tors, which we choose to write in the form (equivalence with the list in |35j 
can be proved by using Fierz rearrangement) 



(A.9) 



In closing this section it is important to note that the form of the Symanzik 
effective Lagrangian enjoys some interesting degree of universality in the 
sense that formally its zero mass limit (m — M c e r = ji q = 0) only depends on 
discretization details, like the form of the gauge action, the specific expression 
of the lattice derivatives, or the value of csw> but not on whether one shall 
be eventually dealing with standard or twisted Wilson fermions. 

When quark mass terms are switched on, the QCD vacuum gets polar- 
ized driving spontaneous chiral symmetry breaking. The information about 
spontaneous chiral symmetry breaking effects is embodied in the Symanzik 
analysis by assigning to the continuum correlators in the expansion appropri- 
ate values consistent with the residual exact symmetries of the target QCD 
theory. In this way the relative "orientation" of the quark mass term with 
respect to the explicitly chirally breaking terms of the Symanzik low energy 
Lagrangian (describing the effects of the Wilson term present in the lattice 
action) becomes crucial for determining the size and the order in a of cutoff 
artifacts in correlation functions (TTJ H] . 



A. 2 Critical mass in Wilson fermion LQCD 

Both in the case of standard (fi q = 0) and twisted mass (fi g = 0(a°)) Wilson 
fermions the critical mass is taken as the value of mo at which the PCAC mass 
vanishes. It is, however, useful to examine separately the two cases, since 
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the resulting structure of lattice artifacts affecting these two determinations 
of the critical mass will be significantly different, as alluded to at the end of 
the previous section. 



A. 2.1 Critical mass in twisted-mass LQCD 

In twisted-mass LQCD the condition (which we write in the physical ip- 
basis (EOl) ) 

J2(^loT 2 ^)(x,t)^ l5 r^)(0)) = (A.10) 



a 



leads to a determination of the critical mass which is "optimal" (M° r pt ) in the 
sense that with this choice all the leading chirally enhanced cutoff effects (i.e. 
those of relative order (a/[i q ) 2k , k = 1,2, .. . with respect to the dominant 
term in the a — > limit) are eliminated from lattice correlators [I]. 

In the spirit of the Symanzik expansion the condition (lA.lOp must be 
viewed as a relation holding true parametrically for generic values of a (and 
fi q ). As a consequence it is equivalent to an infinite set of equations, where 
each equation results from the vanishing of the coefficient of the term pro- 
portional to a k , k = 0, 1, 2, . . ., in the expansion of the l.h.s. of fTA~T0]) B 

That the condition flA.101) yields an acceptable estimate of the critical 
mass can be inferred from the observation that the first of these equations 
(the one which corresponds to the vanishing of the a term) is nothing but 
the continuum relation 

rf 3 x((^ 7 orV)(^t)(^75rV)(0)) =0, (A.ll) 

cont 

by which restoration of parity and isospin is enforced in the lattice theory. 
This means that, if (the 0(a°) term in) tuq is chosen so as to fulfill eq. (lA.lOp . 
then we will simultaneously have m = in eq. ( 1A.5I) and the identification 
on the lattice of the operator -?/>7ot 2 ?/> with the time component of the vector 
current, Vq. At this point to more clearly understand the further implications 
of eq. flA.lOp . it is convenient to explicitly write down the first few terms of 
its Symanzik expansion for which one gets 

= - 



£x J d 4 y(Cf tm (y)V 2 (x)P 1 (0)) 
+ a [ d 3 x([A 1 V 2 (x)P 1 (0) + V 2 (x)A 1 P l (0)]) 



9 The practical problems in determining mo from (|A.10|) . which are related to subtleties 
associated to the exchange of continuum (a — > 0) and chiral (fi q —> 0) limit and to the 
possible limitations of the validity of the Symanzik expansion of the correlator (|A.10|) . 
have been discussed at length in refs. [H [3S] and will not be repeated here. 
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a 2 d'x <fy(Cf^{y)V Q \x)P\Q)) 



cont 



+ a 2 |d 3 a;([A 2 \/ 2 (x)P 1 (0) + \/ 2 (x)A 2 P 1 (0)])| cont + O(a 3 ). (A.12) 

Eq. ( 1A. 12h has been written under the assumption that the 0(a°) piece of mo 
has been chosen so that condition flA.110 is fulfilled. As a consequence (see 
eq. (1A.7P ) £j^ tm is a parity-odd and flavour breaking operator of the form 



/'Mtm 
*"5 



h-sw^l^a ■ Ftp - S^^ipi^r 3 ^ + 0(ji q ) , (A.13) 
while £^ tm (see eq. (IA.8[) ) can be split into the sum of two contributions 



/'Mtm 
^6 



P— even 



^Aqcd^Tst 3 ^ + 0{ii q ) 



(A.14) 



with £^- even a parity -even operator. With the notation AjO in eq. (1A. 12|) 
we indicate the operators of dimension do + j that correct to 0(a 3 ) the 
(local) operator O. Their parity is equal to (— 1) J times the parity of O (see 
Appendix in [1]). We explicitly remark that in eq. ( 1A.12I) we have omitted 
the 0(a 2 ) terms with two integrated insertions of C^ tm and the associated 
contact terms, because they all vanish by continuum parity. Finally for the 
purpose of the present argument we have ignored operators having coefficients 
proportional to fi q . 

At 0(a) the condition implied by eq. ( 1A.12I) reduces to [4 J 







d 3 x / dSdh.sw^l^a ■ FiP 



5 1 A 2 QC ^t l5 r^}(y)V 2 (x)P 1 (0)) +0(/i 



cont 



"qj ' 



(A.15) 



since the first term in the second line of eq. ( 1A.12I) gives a vanishing con- 
tribution to eq. (IA.15P in the chiral limit and there are no 0(a) operator 
corrections to P 1 (i.e. AiP 1 = 0). 

In the chiral regime, where intermediate pion states dominate, one recog- 
nizes that eq. ( 1A.15I) leads to the condition which determines the 0(a) piece 
of the optimal critical mass. We recall that with the definition 



^ = ^(^) = (^lE« 2 ^ + m 5 k 3 (o)) 



cont 



which to leading order in a reads ^ ~ tm 1 7r 3 (0) ) | con t, eq 

comes, in the chiral regime fi q <C Aq CD , the constraint [4] 



& = > 3 (0)> 



cont 



+ 0(a 2 



0(fj, q 



0(a 2 



(A.16) 
HE) be- 

(A.17) 
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which again should be looked at as an infinite set of equations, one for each 
power of a, with the leading one fixing 5\ in terms of b 5] sw- 

At 0(a 2 ) since by continuum parity invariance one can derive the equation 

J d 3 x J d 4 y(£ p 6 - c ™(y)V 2 (x)P\0))\ cont = 0, (A.18) 
J d 3 x([A 2 V 2 (x)P 1 (0) + V \x)A 2 P\0)])\ cont = , (A.19) 

while the terms with the insertion of i/ji^r 3 ^ do not vanish 

'd 3 x J A(^75r 3 ^(i/)K ) 2 (x)P 1 (0))| cont ^ 0, (A.20) 

the condition implied by eq. ( 1A.10I) (or eq. ( 1A.12I) ) yields 5 2 = in ( 1A.14I) . 
The important consequence of this argument is that the estimate of the 
critical mass provided by ( 1A.10I) is not affected by 0(a 2 ) corrections. 

The above line of reasoning can be generalized to all orders in a lead- 
ing to the conclusion that, if the critical mass is determined in twisted- 
mass LQCD by means of the condition (lA.lOp (whose Symanzik expansion 
takes the form ( 1A.12I) ). the expansion of the critical mass will only display 
0(a 2p+1 ), p = 0, 1, . . ., lattice corrections, with coefficients 5 2p+ i determined 
by constraints analogous to (IA. 15|) . 



A.2.2 Critical mass in standard Wilson fermion LQCD 

In the case of standard (clover) Wilson fermions the condition for the van- 
ishing of the PCAC mass is {b = 1, 2, 3, no sum over b) 



2E,(P b (x,t)P b (0)} 



"iPCAC 



for /i q = 0, (A.21) 



which is analogous to eq. ( 1A.10I) with the crucial difference that now [i q = 0. 
The condition ( 1A.21I) is in practice implemented by looking for the limiting 
value of m at which mpcAc ~~ > + . We will call the estimate of the 

critical mass obtained in this way 0. 

Setting mo = M^ r w in the standard Wilson action corresponds to have 
the continuum mass in £4 equal to zero. At the same time lattice artifacts of 
order a k , k = 2,3, ... in will be correspondingly described by terms of 
the kind a k §l w Aqq^xx in L^. The 0(a) term is absent, if the lattice theory 



10 This non-perturbative estimate of the critical mass is the one that is implicitly adopted 
in the studies of standard (clover) Wilson LQCD where the renormalized quark mass is 
defined as m = Z y i^p 1 mpcAC- 
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is clover improved. For the rest, unlike what we have shown to happen in 
twisted- mass LQCD (see sect. I A . 2 . 1 j) . discretization errors of any order in a 



will affect the critical mass determination provided by eq. ( 1A.21I) . 



As in the case of twisted-mass LQCD, also here one must look at the 
conditions coming from the vanishing of the PCAC mass as equations that 
fix the values of the coefficients 8^ w in terms of the matrix elements of Ck 
between one-pion states. For instance, at 0(a 2 ) from the symmetries of the 
standard Wilson action (and the form of the associated Symanzik expansion) 
one cannot conclude anymore that vanishes. Rather the value of this 
parameter is fixed through eq. (1A.21I) by the condition 

(7r(0)|£ 6 |7r(0)) =0, (A.22) 

cont 

where Cq is the full six-dimensional operator of the Symanzik effective La- 
grangian (given in eq. ( 1A.8I) ) which, we recall, also includes the (^A-qcdXX 
term. Eq. (IA.22P follows from the fact that the lattice pion squared mass 
and mpcAcU by construction vanish at the same value of rriQ. Furthermore 
these two quantities are found to be linearly proportional to a very good ap- 
proximation in the region where they are both small, implying that possible 
additive cutoff artifacts affecting mpcAcli and the squared pion mass are 
also proportional to each other. 



Appendix B - 0(a 2 ) corrections to the squared pion 
mass 



The correlator of interest for extracting the squared pion mass is the four- 
dimensional Fourier transform of the two-point subtracted pseudoscalar 
correlator (no sum over b = ±,3) which reads 



(B.l) 



T h L {p) has a pole at p 2 
eq. (EH])) 



2 

-m b 



l with a residue given by |Gv,||, where (see 



(n\P b \n\0)) 



L 



(B.2) 



11 Since we are interested in determining the structure of the cutoff corrections affecting 
the lattice pion mass, we can always imagine that contact terms, which do not display 
the pion pole, have been subtracted out. In the formulae that follow the superscript 
usubtractedn j g ^ nug a i wa y S understood. 
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To proceed further we need to write down the Symanzik expansion of ^l(p) 
up to order a 2 included. This gives 



Y,e ipx (P b (x)P b (0))\ L = J d 4 xe ipx (P b (x)P b (0)) 
d 4 xe ipx (A 1 P b (x)A 1 P b (0)) 



cont 



+ 
+ 

+ 



d 4 xe zpx {A 1 P b (x)P b (0) I d 4 yCf tm (y)) 



cont 

4„, rMtrni 



cont 



d*xe ipx (P b (x)A 1 P b (0) J d 4 yCf tm (y)} 
d 4 xe ipx (P b (x)P b (0) J d 4 yCf tm (y) J d 4 y' C 



Mtm/ 
5 



— a 



d 4 xe ipx (A 2 P b (x)P b (0)) 



cont 



— a 



d 4 xe ipx {P\x)P\<S) / A 4 (2/)) 



M> cont 

d 4 xe ipx (P b {x)A 2 P b {0)) 
0(a 4 



cont 



(B.3) 



The absence of odd powers of the lattice spacing in this formula is guaran- 
teed by the results of ref. [TT], as we are working at maximal twist. We will 
also assume that the critical mass has been determined in the optimal way 
described in sect. IA.2.11 so that by setting mo = M° r pt the condition (1A.17jl 
holds true. We also recall that the operators A\P b and A 2 P b are the di- 
mension four and five terms that are needed to compensate for the 0(a) and 
0(a 2 ) operator corrections arising from the contact of C^ tm (y) and C^ tia (y), 
respectively, with the operators localized at the points x and 0. 

Since we are interested in the lattice corrections to the pion mass, we must 
look in the r.h.s. of eq. ( IB. 31) for the terms which may display the double pion 
pole factor (p 2 + m\Y 2 . 

Among the many terms in the r.h.s. of eq. (IB.3[) only the fifth and the 
last display the double pole we are after. Putting together eqs. (12.111) . (IB.2|) 
and ( IB. 31) and taking the limit p — > at small m 2 (but always parametrically 
larger than a 2 [1]), from the ensuing pion pole dominance we arrive at the 
formula 



\G w s\ 



2 

m 3 



\G n \ 2 




L ml 


cont . 



I -a 2 



^ 3 (0)|^5 5 (0)k 3 (0)> 



mi 



cont 



O(-t), (B.4) 



mi 



with 



/'Mtm 
*-6Sz55 



(0) 



£f m (0) ~\j d 4 xCf tm (x)C 



Mtm 
5 



(0). 



(B.5) 



12 For short from now on we will simply write m 2 (with no isospin index as in the 
continuum target theory all three pions are mass-degenerate) for mJ| cont . 
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A simple Taylor series resummation leads precisely to eq. (12.1 1) of the text, 
which we stress also agrees with the results derived in x?T [T21 E21 EH 38J . 

We close the Appendix by noting that the absence of a contribution from 
(7r ± (0)|£|^ tm (0)|7r ± (0))| C ont in eq. (12.61) is a consequence of the commutation 
relation o 

from which one gets the SPT 



^(0)1^(0)1^(0)) 



COIlt 



-I 



(Q\[Q% £^(0)11^(0)) _ + OK) = OK) . (B.7) 



COIlt 



Appendix C - Expression and estimates of A^ 5 ' 3 

In this (long) Appendix we present arguments showing that Ag 5 , b — 3, ± 
(eq. (I2.5p ) yield contributions to rn^ b \i, (eqs. (I2.6P and (12. 7p ) which are para- 
metrically (in the limit m 2 — > 0) and/or numerically negligible, provided 
Mtm-LQCD is defined by setting m (see eq. ( 1A.1I) ) at its optimal critical 
value determined by the condition (1A.10I) . 

In sect. IC.ll we provide the expression of Ajj 5 and A^ 5 (eq. (12. 5p ) based on 
standard field theoretical manipulations, i.e. insertion of intermediate states, 
or alternatively LSZ representation and reduction formulae, and, wherever 
applicable, SPT's [281 EH]. In sect. IC.2I we give a numerical estimate of 
the magnitude of A| 5 and A^ 5 for the typical lattice volume and pion mass 
values relevant for the recent ETMC simulations. The details of our estimate 
of the key matrix element £ w (eqs. (IA.16I) and (1A. 17|) ) . which controls the 



numerically dominant (in the small m 2 region) contributions to A55 , are 
deferred to sect. IC.31 



C.l The theoretical analysis of A| 5 and A^ 5 

For the purposes of our study we restrict attention to Nf = 2 QCD with 
light u and d quarks yielding light, weakly interacting pions. In this regime, 
on general grounds, one can write in the infinite volume limit the spectral 
decomposition of the continuum quantity A| 5 (eq. (12. 5p ) in the convenient 

13 In Mtm-LQCD the flavour symmetry group is the "oblique" SU(2) D b group generated 
by the charges Q\, Q~ A and Qy. These conserved charged are associated with the (one- 
point split) currents reported in sect. 4.1 of ref. [H] . In the notation of that paper we are 
dealing here with the case r = 1, u> r = ir/2. 
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form 1^1 



A 55 = -^(7T 3 (0)|£f m K 3 (0)7T 3 (0)) + (C.1) 



mi 



^55 - - 2^ m 2 _ m 2 + ■ ■ • ' y L - 2 > 

where we have explicitly separated out from the remainder, R^ 5 , the con- 
tribution of the semi-disconnected one-pion pole (sometimes referred to as 
"tadpole" in the following) and the contribution of the cut over two-pion 
states below the inelastic threshold. In order to have more manageable for- 
mulae we have decided to reduce the initial integrations over the momenta of 
the two pions to the integration to the single center of mass energy variable, 
and we have set 



k{E) = V(£/2) 2 - ml . (C.3) 

The remainder in eq. ( 1C.2I) comprises a sum over one-particle scalar states 
with zero isospin 0, as well as contributions from heavier and/or more com- 
plicated multiparticle states which we have simply indicated by ". . ." . Our 
conventions are such that in eq. OCIO the single-particle meson states a 
(a = 7r, <7°, . . .) are introduced with the Lorentz covariant normalization 
(b,b' = 1,2,3 are SU(2) isospin labels) 



(a\p)\a b '(p>)} = 2E a (p) (2^)%^ - f) , E a (p) = y/p + m*. (C.4) 

The states |"P^ 3=0 ^7T7r, E) are two-pion states with total energy E, zero total 
three-momentum and zero total isospin third component. In particular the 
symbol V^ =0 ^ denotes the isospin projector onto the two-pion states with 
vanishing total third component (1% = 0). Here we recall that two possible 
states contribute which one has to sum over, namely the state with two 
neutral pions and that with one negatively and one positively charged pion. 

Similarly the expression of the (infinite volume limit) spectral decompo- 
sition of A 5 5 (eq. ( 12.51) ) will read 

A± = -^(^(0)1^-1^(0)^(0)) + (C.5) 



14 For short in this Appendix we will drop the subscript | con t from continuum QCD 
quantities. Lattice quantities will still be labeled by the subscript \l. 

15 Phenomenologically the lightest of these states should be identified with the ao(980) 
meson. It should be kept in mind that, since in the actual simulation data we will be 
considering the lightest "pion" mass is about 300 MeV, the lattice state corresponding to 
the "ao (980)" -particle is expected to have a mass somewhat above 1 GeV. 
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167T 3 



^55 = - 2 _ w 2 + • • • ' I " 6 ) 

with P^- 3) the isospin projector onto the state of two pions having third 
component ± and 3. 

On general grounds one expects the numerically dominating contributions 
to to come from the interaction of the propagating pion (the initial and 
final particle in eq. (12.51) ) with a neutral pion created from the vacuum, i.e. 
from first term in the r.h.s. of eqs. (IC.ip and (10. 5[) . In the next sections we 
will show that this is indeed the case. The important observation is that in 
Mtm-LQCD with optimal critical mass these dominant tadpole contributions 
are 0{m\) and moreover equal for A 3 5 and A 55 up to (very small) O(m^) 
terms, hence they get canceled in the pion squared mass difference. The 
estimated leftover correction coming from the elastic two-pion cut as well as 
the remainders (eqs. ( 1C.2I) and ( 1C.6I) ) will be shown to be tiny. 

We now examine in turn the various terms contributing to A 3 5 and A^ 5 
and discuss their parametric behaviour in the small-m^ regime. We shall see 
in particular that for symmetry reasons the contributions from the various 
intermediate states to Af 5 tend to be negligible compared to those contribut- 
ing to A 3 5 , with the exception of the tadpole term and the one coming from 
the lowest-lying two-pion intermediate state. 

C.l.l Tadpole 

That the first terms in eqs. ( 1C.1I) and ( 1C.5I) are equal to leading order in the 
chiral expansion follows, e.g. by reducing the neutral pion in the ket state, 
yielding | 
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(7r 3 (0)|£f m |7r 3 (0)7r 3 (0)> = (vr^O) |£f m | 7 r 3 (0)7r ± (0)) . (C.7) 

An explicit estimate of the tadpole terms can be given on the basis of the 
SPT relations 

(^(0)1^1^(0)^(0)) S ^ T --^(0|£f > 3 (0)> = -| , (C.8) 

J 7T J IT 



from which one gets 



A 3 

^55 



tad rn? 



%< 7 r 3 (0)|£f > 3 (0> 3 (0)> S = 
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In this Appendix with the symbol ~ we denote equality up to 0(a 2 ) corrections, with 

SPT SPT 

= equality up to 0(m 2 ) terms and with ~ equality up to 0(a 2 ) and 0(m 2 r ) corrections. 
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A^5 



A = — %<7r ± (0)|/:^ t m |7r 3 (0)7T ± (0)) S ^ T . (C.9) 

tad m 2 N 5 /2 m 2 / 



We remark that when the untwisted mass, mo, is fixed at its optimal value, 
eq. ( 1A.17I) holds |4J, implying 

- -oK), (c.io) 



while in the opposite case one finds that this quantity tends to increase as 
we go towards the chiral limit, leaving large 0(a 2 /m 2 r ) cutoff effects in the 
neutral as well as in the charged lattice pion mass r 7 l . 

C.1.2 Two-pion cut below the inelastic threshold 

In a finite three-dimensional volume, V = L 3 , the integral over the two- 
pion cut can be rewritten as a sum over discrete states through the formal 
replacement 

^ k(E) 1 k(E e ) 

2m, E*-Tr$'" f^ Pv {E t )Ej-ml---' ^ ) 

where £ max is the number of allowed two-pion levels that can be hosted in 
the volume V below the inelastic threshold (4m„-). The relation between the 
allowed values of i and the level energies is established combining eq. (1C.3I) 
with the formulae [391 HO] 

kl 

£tt - 6(h) = 4>{q) , q = —, £>0, (C.12) 

where S(k) denotes the s-wave 7T7T phase-shift in the appropriate isospin chan- 
nel and <p(q) is a known j3U] kinematical function. Finally (see e.g. ref. 



pyy E ) ^ E > (O- 13 ) 



is the two-particle state density with 8'{k) = dd(k)/dk and 4>'(q) = 
At the volumes (L = 2.1 and L = 2.8 fm) that will be considered in sect. IC.2I 
and for the lightest pion mass (m n ± ~ 300 MeV) one has / max = 2. 



17 This result is in agreement with the finding of theoretical studies [42l[T3l|4] and indi- 
rectly confirmed by the observation of large, positive 0(a 2 ) lattice artifacts in the quenched 
Mtm-LQCD computations [121 [H] of the charged pion mass and the corresponding decay 
constant (computed from the WTI relation / ff ± = 2jj, q G^/nn? r± ). 
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• I = - The contribution of the i = term in the sum fjO. 11 j) can 
be computed in an almost model-independent way observing that the lowest 
energy level corresponds to have the two pions at rest with Eq = 2m n 0. 
Making use of the SPT's (1C.71) - (1C.8P and taking into account the presence of 
the isospin projectors V^ = ^ and V^'^ in eqs. ( 1C.1I) and ( 1C.5I) . respectively, 
one finds 



A 3 

^55 



3« 



1=0 



[El - ml) ft 



{El -ml) ^ 



i?o=2m, 7r 



i?o=2m, 7r 



where 



hi 



k(E ) 



8n 2 p v (E ) 



(l + 0(l/L) 



(C.14) 
(C.15) 

(C.16) 



E =2m n 2m n V 

will be in the following approximated by its infinite volume limiting value, 
(2m 7T V)~ 1 , which we remark is independent of the expression of the phase 
shift S(k) appearing in p v (E). 

Two observations are in order here. The first is that the quantities (10. 14[) 
and ( 1C.15I) only differ because different isospin combinations of two-pion 
intermediate states contribute. The second is that, as in the tadpole case 
discussed above, both contributions are of O(m^), if eq. flC.101) holds, i.e. if 
the critical mass has been set to its optimal value. If this is not so, large and 
unequal, 0(a 2 /m^) cutoff effects will be left out in the neutral as well as in 
the charged lattice pion mass. 

• i > 1 - As we said, given the magnitude of L we consider, in the energy 
interval 2m n < E < Am n two-pion states with I = 1 and 1 = 2 will contribute 
with terms of the form 



^55 



2h\ 



mi 



{7r b {6)\Cf tm \V. 



(...) 



7T7T, En) 



k{E e 



87T 2 Pv (E e 



(C.17) 



with b equal to 3 or ± and standing for V^ = ^ or V^' 3 \ respectively. 
Evaluation of A| 5 |£ = i 2 is possible only at the price of making some reason- 
able assumption about the energy behaviour of the pion matrix elements in 
eq. (lC.17p . The simplest thing is to assume that these matrix elements are 
constant in energy in the range 2m n < E < Am^. The quantities h\ in 
eq. ( 1C.17I) can be straightforwardly evaluated in terms of k{E^) and py{Ei) 
for the I- values of interest, once a reasonable parameterization of the phase- 
shift 5{k) (see sect. IC.2j) is inserted in eq. f)C.13l) . 



18 Actually there is a finite size correction proportional to the (J = 0, / = 0) 7T7T 
scattering length [35] ■ To be precise one gets Eg = 2171^ — 47ra[} (m^V)^ 1 (1 + 0(1/L)). 
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Concerning Af 5 |^ =1 2, the important remark is that these quantities in the 
chiral regime are suppressed by a factor of m\ with respect to their neutral 
counterparts Ajy^i^. The result follows from eq. ( 1C.17I) by noting the SPT 
relation (actually valid for all t > 0) 

iU (7^(0) \£f* >*\Vg*>mc,E e ) = (n\[Qi, £f m ] \V^hn,E £ ) = 0(ml) . (C.18) 

As a consequence, assuming the validity of this SPT for m n ~ 300-r-400 MeV, 
we shall neglect Af 5 \£ =lj2 with respect to Ajy^i^ in our subsequent numer- 
ical estimates. 

We also observe that, if £j^ tm itself is an operator of 0(/x 9 ) = O(m^), as 
it happens in the LO-^PT description of Mtm-LQCD with optimal critical 
massif, the r.h.s. of eq. ( 1C.18I) is actually an O(m^) quantity for states with 
E e = 0(m T ), while (n 3 (Q)\£™ tra \V£ =0) 7rrr, E e ) is 0(m£). In this situation 
one thus finds A S5 |^ =1) 2 = 0(m®) and A| 5 |£ =12 = O(m^). The resulting m\ 
suppression factor is in agreement with the general m^-behaviour (implied by 
the SPT (IC.lSp ) of the charged pion contribution with respect to the neutral 
one. 

C.1.3 The remainder 

The remainder comprises the two-pion cut integral from 4m n to oo and the 
contribution from all other possible single and multiparticle intermediate 
states. 

The first important observation is that all the contributions to Rf 5 are 
of O(m^). This result can be derived in close analogy to eq. (1C.18I) above, 
namely by reducing the charged pion in the matrix elements (see eq. ( 1C.6I) ) 
(7r ± (0)|£| Itm |a ± ), where |a ) is a state with a non-zero energy in the chiral 
limit (e.g. the one-particle scalar state appearing in eq. (10. 6[) ) . One gets, 
in fact 

^<^(o)|/;f n V> = (n\[Q%£T m ]\o^) = oK). (c.19) 

A similar result does not hold for R^ 5 because the commutator [Q% C^ tra ] 
does not vanish in the chiral limit. Based on eq. ( 1C.19I) we shall neglect the 
term Rf 5 with respect to i?f 5 in the numerical estimates below. 

19 In LO-xPT there is only one operator (see e.g. ref. |13j ) with the same quantum num- 
bers as Cf tm (eq. (|Xl3|) ). Hence the condition £ w = (tt\£f tm \ir a (0)) = 0(fi g ) implies 
that the unique effective operator representing at LO of xPT must itself be propor- 

tional to the quark mass. Obviously for matrix elements involving states with increasing 
rest-frame energy E the LO-xPT description progressively looses its validity. 
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In this context it should also be noted that the sum Af 5 \ i= i+Af 5 \ e=2 +Rf 5 , 
being a negative quantity, yields a (tiny) positive contribution to the phe- 
nomenologically negative (2j [3j squared pion mass splitting. This term (see 
eqs. (I2.6p - (12.7p ). if included in A^ 5 , would go in the direction of increasing 
the negative contribution of (, n that is necessary in order to reproduce the 
observed value of 771^3 1 2, — In- 
coming back to the evaluation of the remainder terms, we shall hence 
limit below our discussion to i?| 5 . The latter, we recall, includes, besides 
terms from one-particle scalar intermediate states explicitly shown in eq. ( 1C.2I) . 
also a cut contribution of the form 

as well as negligible contributions from more complicated (and heavier) mul- 
tiparticle intermediate states. The cut contribution above can be estimated 
by assuming the standard Kallen-Lehmann (large energy) 1/E 2 behaviour 
for the modulus square of the C^ tm matrix elements and by setting the coef- 
ficient factor in front to a value that matches the contribution from the £ = 2 
two-pion intermediate state. As for the one-particle scalar state contribu- 
tion, we make the plausible assumption that it is of the same (very small, see 
sect. IC.2p size as the cut contribution, and simply double the latter in order 
to get our numerical estimate of i£f 5 . 

Terms from multiparticle intermediate states of increasing mass are ex- 
pected to give tiny and progressively smaller and smaller contributions, owing 
to explicit energy denominators and the asymptotic high-energy behaviour 
of the matrix elements. 



C.2 Numerical considerations 

In this subsection we want to provide a numerical estimate of the order of the 
magnitude of the difference A^ 5 — Af 5 relying on the results of the previous 
subsection and the numerical evaluation of £ n given in sect. IC.31 

We shall see that A^ 5 , which according to the considerations of sect. IC.ll 
is an O(m^) quantity, takes a value much smaller than a _2 Am^.| L ltm . Numer- 
ically Ag 5 is of similar siZG ctS ctS the parametrically 0(m°) contributions 
from i?5 5 are tiny and actually smaller than the O(m^) terms. 

These facts imply (see eqs. (12. 6 p and (12.71) ) that lattice corrections to the 
squared charged pion mass are also small and of order a 2 m 2 (or a 4 ), while the 
main contribution to the lattice artifacts of the squared neutral pion mass 
comes from as announced in the text (sect. [2D- As we already said, in the 
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difference A| 5 — Af 5 the contribution coming from the tadpole terms exactly 
cancels. 

Using the value ( 1C.54I) of £ n and the estimate h% = (2m n V)~ 1 , we can 
provide an explicit numerical evaluation of the three contributions to A 3 5 
and A^ 5 described in sects. IHXT1 EU and IHX31 To arrive at this result 
we also need to know the energies of the few lowest tttt levels living in our 
finite simulation box at the actual value of the simulated pion masses. In 
the situation corresponding to the lattice data used for our best estimate 
of £„-, we have m n ~ 300 MeV and L ~ 2.1 fm. We checked that results 
do not change significantly with the volume by repeating the same analysis 



for L ~ 2.8 fm. In order to compute kt and hence Eg = 2ykf + from 
eq. (1C.12I) . as well as pv{Ee) from eq. ( 1C.13I) . we parameterize the s-wave im 
phase shift in the form (suggested by LO-^PT |45j ) 

with f = 86 MeV being the chiral limit value of the pion decay constant 
(as obtained from recent ETMC studies [31 [2])- The formula (lC.20p (with 
m n at its physical value 140 MeV) describes rather well the experimental 
7T7r-scattering data for center-of-mass energy E in the range 2m n to Am^ (see 
e.g. the phase shift data compilation in ref. [43]). 

Following the discussion in sects. IC.1.1[ 10. 1 .21 and IC.1.3[ in terms of the 
reference scale 03=3.9, whose value is currently estimated [3] to be about 
2.3 GeV, we obtain the following numerical results. 

• - Tadpole 



A 3 

^55 



= ^55 

tad b0 



SPT 



2 

+0.00032 



~ 0.00050_™a^ 3 .9, (C21) 

taC1 JTV m TV 



where the quoted error reflects the uncertainty on £ n from eq. 00.540 . 

• - Two-pion cut below the inelastic threshold. For the I = term one 
obtains 



A 3 

^55 



A^ 5 



= "HI ~ -(°- 000063 -o:o=) «A. 9 , (C22) 

~ -(°- 000042 -o:=) hU, • (c.23) 

7T J 7T 



For the i > levels one gets smaller and smaller contributions. A reasonable 
estimate of their order of magnitude is — 0.000020al =3 9 and — 0.000015a^ =3 9 
for A 3 5 |^ at £ = 1 and £ = 2, respectively, to which we attach a generous 



26 



50-60% relative uncertainty. We neglect here A^ 5 |^ =12 , which according to 
SPT arguments, are significantly smaller {m\ suppressed) than their neutral 
counterparts. 

• - Remainder. According to the results of subsect. IC.1.31 we can limit 
consideration to K\ % . Despite the difficulty of evaluating all the terms con- 
tributing to i?| 5 , following the procedure discussed in sect. IC.1.31 we can 
reasonably estimate it to be about — 0.000080a^l 3 9 , with, as before, a 50- 
60% relative statistical uncertainty. We thus find that R\ h is of the same 
order of magnitude as Aiy^o or Y^=\ A55 \e_. 

• - Total. Putting everything together we finally get 

Aj? 5 ~ +(0.00032 ± 0.00019 ± 0.00008)a^ 3 9 , (C.24) 

A| ~ +(0.00046 ± 0.00027 ± 0.00005)a^ 3 9 , (C.25) 

and 

A^ 5 - A| ~ -(0.00014 ± 0.00008 ± 0. 00004) a^l 3 .9 • (C-26) 

The two "errors" quoted in the equations above reflect the uncertainty (dom- 
inantly of statistical nature) on the value of £„- from sect. IC.3I and the sys- 
tematic indetermination inherent to the assumptions we made above about 
the energy behaviour of the matrix elements of £5 m , respectively. 

For the sake of comparison here we only recall that (as one can infer 
from Fig. [Hand the value ro/al/3^.9 ~ 5.2) at ~ 300 MeV direct lattice 
measurements yield a| =3 9 Am^|^™ 3 9 ~ —0.0067(20), i.e. a number much 
bigger than ai =3 9 Ag 5 , or a^ =3g Af 5 . Comparisons of this kind are performed 
and discussed more thoroughly in sect. I2.2L 

C.3 Estimating ^ from lattice two-point correlators 

We describe in this section the analysis, based on the Symanzik expansion, by 
which we arrive at a numerical estimate of the crucial continuum quantity ^ 
(eqs. ( 1A.16I) -( TA.17I) ) that has been employed in sect. IC.ll to parametrize the 
various terms contributing to A^ 5 and A| 5 . Since in this numerical analysis 
we will be exploiting Nf = 2 Mtm-LQCD data at sufficiently small m n ~ 
m^± \l around 300 MeV, we feel safe to use SPT's to relate among themselves 
some of the (continuum QCD) matrix elements that appear in the Symanzik 
description of lattice quantities. The SPT approximation is correct up to 
O(m^) relative corrections, whose relevance is checked by going to heavier 
quark masses up to m n ± \l < 450 MeV. By extending the analysis from data 
at lattice resolution a" 1 = 2.3 GeV {(3 = 3.90) to those at a -1 = 2.9 GeV 
(j3 = 4.05) and comparing (at the lowest-lying pion mass) results at box 
linear size L ~ 2.1 fm with those at L ~ 2.8 fm, we have checked that within 
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our (large) statistical errors the numerical estimates relevant for this paper 
display no significant discretization and/or finite size effects. 
The key correlator we consider is P 



C L PS (t) = £~<P 3 (f, t)S\0))\ L - L 3l -(P 3 (0))(S°(0)) 



t>0, (C.27) 



which in the continuum limit is an 0(a) quantity because it violates par- 
ity (and isospin) invariance. The Symanzik expansion of its renormalized 
counterpart reads 

Z s0 Z P C L PS {t) = a(J dxP 3 (x,t) J d 4 yCf tm (y)±S°(0)) + 0(a 3 ) = 



m 2 O 
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~2 



-m w t 



2m„ 



+ 



aG n (n 3 \£f tm \a e{i )(a eS \i-S»\n)- 



ni 



°"cff 



mi 



2m n 2m„ 



+a\(7r 3 \P 3 \a cS )(a cS \l-S \n}' 
mi 2 



2m, 



+ 0(a 3 ) 



+ 

(C.28) 



°"cff 



where P 3 and S° denote the renormalized operators P 3 and S°, respectively, 
and G n was defined in eq. ( 12. lip . The symbol |cr e ff) represents the lightest 
(isosinglet) scalar state that can be created from the vacuum by the operator 
S°. For the quark mass values and in the t-range of interest for the present 
analysis of Cp S (t) (because of the unfavourable signal-to- noise ratio inherent 
in the evaluation of the quark disconnected piece of the correlator), we are 
in fact sensitive within errors only to the neutral pion and \cr e s) intermediate 
states. In a finite box and for sufficiently small quark masses, the symbol 
|cr e ff)(cr e ff | in the r.h.s. of eq. (1C.28h is to be interpreted as 



|cr e ff)(oeff| = hlY^ IAV^I^ttV^I 



m a cti — 



2m„ 



(C.29) 



6=1 



In the following we shall write simply \a) for |cr e ff). 

The three terms in the r.h.s. of eq. (1C28h come from the three different 
possible time-orderings of the inserted operators. We have only retained the 
pole contributions from the two lightest states and the associated partially 
disconnected terms that must go along with them in order to comply with 
Lorentz invariance. Since in Mtm-LQCD with optimal critical mass all the 



20 We recall that all the correlators we shall be dealing with are "connected" in the strict 
sense of Quantum Field Theory, though some contribution only through gluon exchanges. 
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leading terms (linear in a) in eq. (IC.28I) are of the same order in m n , we have 
ignored here the 0(a) lattice corrections to the operators S° and P 3 , because 
they would give contributions to the Symanzik expansion suppressed by an 
extra (relative) factor m\ and are numerically small for the considered 
range of quark mass and t-values. 

With the help of eqs. flCT8]) and (jC^Sli as well as the SPT relations 

' " ^V) 8 ^!^, (C30) 



(n 3 \P 3 W) S =-yjn\-S»\a), (C.31) 
the Symanzik expansion ( 1C.28I) can be rewritten in the simpler form 



Z P Z s oCp S (t) = T a 



2m n f n ml 



3h 2 



m„ 



-m,yt 



mi 



mi 



fn m l 
+ 0(a 3 ) 



mi 



+ 



(C.32) 



2m a f n ml fl 

where, as we shall see below, all the quantities, but £ n and the negligible 
0(a 3 ) terms, can be directly evaluated from lattice data. Solving, in fact, 
the resulting equations yields an estimate of for £„. for each given quark mass 
fig and lattice resolution a. 



C. 3. 1 Evaluating £ n 

Several two-point Green functions in the neutral pseudoscalar (PS) channel 
(where also quark-disconnected diagrams contribute) have been evaluated, 
among which 

C L PP (t) = a 3 E^ 3 (^ 3 (0))| L -^ 3 (P 3 (0))| L (P 3 (0))| L , (C.33) 

X 

Ci s (t) = a 3 ^\(S\x)S\0))\ L - L 3 \(S\0))\ l (S°(0)}\ l (C.34) 



and 



C L PS it) = a 3 ^l(P 3 (x)S°(0))\ L - L 3 1(P 3 (0))\ l (S°(0)) 

X 

C L SP {t) = a 3 Y, \{S\x)P 3 m\ L - L 3 1(S°(0))\ l (P 3 (0)) 



L ' 



(C.35) 
(C.36) 



21 This is seen by noting that AiS° cx a/j, q P 3 , while AiP 3 comprises a term proportional 
to a[j, q S a and another one proportional to F^F^. This last term has 0(fx q ) matrix 
elements between states involving only the vacuum and one or two neutral pions, as one 
can check by means of SPT's. 
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It should be noted that in connected correlators the mixing of the bare op- 
erators P 3 and S° with the identity does not contribute. 

All lattice correlators are simultaneously fitted to simple two-state ansatz 
of the form 



PS 
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exp(— m^it) 



2m, 



+ d 5 di 
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(C.37) 

(C.38) 
(C.39) 



2m n 3\ L 

Comparing the fit ansatz for the correlators Cp S ss PP with the correspond- 
ing Symanzik expansions written to the leading order in a and in the form 
where only the neutral pion and/or the lightest scalar state contributions are 
retained (consistently with the numerical fact that within statistical errors 
only the contribution of these two states is visible), one obtains the following 
identifications 



c 5 ~ (fi|PV(0)> 



di~(n\-S°\a), 



as well as 



rr rv SPT 

ZgoZpCiC^ ~ a 



3h% 



rnz 



fl ml - ml 



(C.40) 



(C.41) 



A similar equation between ZsoZpdid$ and the coefficient of the e~ mat j (2m cr ) 
term in the r.h.s. of eq. ( 1C.32I) can be obtained. This relation is, however, 
of minor interest given the substantial statistical error we have on the fitted 
coefficient d$. Finally the relevant dimensionless ratios Zp/Z s o and 3h^/f^ 
entering eq. (1C.41I) can be estimated via the formulae 



c 5 



Z s o G n ±\ L 



hi spt 1 / di 



(C.42) 



where 



G T± \ L = (n\P ± \n ± (0))\ 1 
P ± = (P 1 ±iP 2 )/V2. 



(C.43) 
(C.44) 



The second of the eqs. (10.421) is a consequence of the SPT (I0.3Q[) and the 
relation 



Zpdi S ~ T v / 3/io( 7I " v 



-S°\ti 3 (6)) 



(C.45) 
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which in turn follows from the identification of a as a two pion state (see 
eq. (IC.29I) ). It should also be recalled that both the renormalized operators 
So and P ± are related to their bare (subtracted) counterparts through the 
renormalization constant Zp, which would then cancel in the ratio di/G^li 
if we were to write it in terms of renormalized quantities. 

Combining together the relations (lC.40p . (lC.4ip and flC42[) . one arrives 
at an explicit formula for £„., namely 



t spt f^mli 



3m; 



3G 2 + in 2 — m 2 - 



-i 



where the charged pion sector 0(a°) quantities m n ±\ L , G n ±\i and 



f n ± = 2/i ? (G 7r ±m 7r 



(C.46) 



(C.47) 



can be evaluated much more accurately (at level of 1% or better) than the 
corresponding neutral sector quantities, m a \ Ll d 1 and C\ (we recall that C\ 
vanishes linearly in a as a — > 0). From the analysis of the neutral pseu- 
doscalar channel correlators specified in eqs. (10.33[) to (10.36j) . we can extract 
m a \L that, within statistical errors (ranging from 10% to 20% depending on 
statistics) turns out to be equal to 2m 7r ±|i at all the values of [i q and at the 
two different lattice resolutions we consider (see Table [2] below) . Using the 
relation flC.47p and setting m a \ L = 2m n ±\ L , eq. (IC.46I) takes the simple form 



or alternatively 



c SPT „ 



c SPT „ 



"iG^± 



h 2 



n -l 



(C.48) 



(C.49) 



The latter expression of £ n is exploited to obtain the numerical estimates 
presented below. 



C.3.2 Numerical results 

We report in Table [2] in lattice units the results for the quantities enter- 
ing eqs. flC.46p - flC.49p . These numbers are obtained from the analysis of 
the correlators flC.33p - flC.36p . as well as their charged pseudoscalar meson 
analogs, evaluated on the gauge configurations taken from the Nf = 2 ETM C 
ensembles with the following bare parameters 



(3 = 3.9 (a 
(3 = 3.9 (a 
P = 4.05 (a 



0.09 fm) , L/a = 24 
0.09 fm) , L/a = 32 
0.07 fm) , L/a = 32 



a/2 q = 0.0040, 0.0064, 0.0085 , (C.50) 
a/i 9 = 0.0040, (C.51) 
afiq = 0.0030, 0.0060 . (C.52) 
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The values of a\i q above are known [31 E] to correspond to quark mass values 
(in the MS scheme at a scale of 2 GeV) in the range from 20 to 40 MeV. 
At (5 — 3.9 and a[i q = 0.0040 data at two different volumes allow to check 
that our estimates of £ n are not significantly affected by finite size effects. 
The results at the two /5-values can be compared by expressing quantities 
in units of r (extrapolated to zero quark mass limit) with the help of the 
estimates [31 [7]) 

r /a\f3 =3 .9 = 5.22(2) , r /a\ p=Am = 6.61(3) . (C.53) 

In order to ease such a comparison, for the case of c±, which is a dimension 
two and O(a) quantity, we also quote the value of cir^a -1 , the latter being by 
construction a dimensionless and O(a ) quantity o. For the key quantity £ n 
we quote for the same reasons r^ n . In Table [2] we also give for each ensemble 
the number of measurements meas.) of the neutral pseudoscalar meson 
channel correlators. 

Inspection of Table [2] shows that the most accurate results are obtained 
for (3 = 3.9, afi q = 0.0040, where the number of measurements of neutral 
pseudoscalar meson channel correlators is the largest and we also have data 
on two different physical volumes. From the results at (3 = 3.9, a\i q = 0.0040 
and using eq. (10.491) . we obtain our best £ n estimate, namely 

& S ~ T -0.00014(4) a^ 3 9 & r 4 ^ S ~ T -0.10(3), (C.54) 

from which eqs. (10.211) (10.26[) follow. One checks in Table [2] that this result 
remains stable within statistical errors at increasing values of the quark mass. 

If we were to use for £ n the expression (10.481) . instead of eq. ( 10.491) . results 
consistent with eq. ( 10.541) would be obtained, though with a larger statistical 
errors due to the fact that di/(y/3G n ±)\L is a quantity of order one with a 
sizeable uncertainty from d\, the effect of which is enhanced in the factor 

(1 — 3G 1 p j . For instance, at /3 — 3.9, afi q = 0.0040 and L/a = 24 (the 

statistically most precise data set) we would get £ w = —0.00015(5) 0^=3.9- 

Results about the various quantities entering eqs. (10.461) (10.491) are also 
shown in Table [21 as a check of their quark mass dependence and scaling 
with a. For instance, in the case of c\ we find by inspection that the values 
of C\r\a~ x display within errors good scaling and no visible quark mass de- 
pendence. This last property is consistent with the use of SPT's we made 
in order to arrive at eq. (10.48|) . The situation for d\ is similar, but with 
somewhat increasing statistical uncertainties at higher quark masses. The 
statistically very precise matrix element G n ±\i,, shows instead a mild quark 

22 The statistical error on cirjja -1 is completely dominated by that on a 2 c\. 
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mass dependence, as its value is seen to change by not more than 10% in the 
quark mass range we have considered. 

Acknowledgments - Illuminating discussions with M. Testa are grate- 
fully acknowledged. We are also indebted to M. Liischer for stimulating 
remarks and P. Weisz for a technical discussion. We wish to thank the rest 
of ETMC for the interest in this work and for a most enjoyable and fruitful 
collaboration. This work has been supported in part by the DFG Sonderfor- 
schungsbereich/Transregio SFB/TR9-03, DFG project JA 674/5-1 and the 
EU Integrated Infrastructure Initiative Hadron Physics (I3HP) under con- 
tract RII3-CT-2004-506078 and by the EU Contract MRTN-CT-2006-035482 
"FLAVIAnet" . We also acknowledge the DEISA Consortium (co-funded by 
EU under the FP6 project 508830) support within the DEISA Extreme Com- 
puting Initiative (www.deisa.org). G.C.R. and R.F. thank MIUR (Italy) for 
partial financial support under the contracts PRIN04 and PRIN06. 

References 

[1] R. Frezzotti, PA. Grassi, S. Sint and P. Weisz, JHEP 0108 (2001) 058; 
R. Frezzotti, S. Sint and P. Weisz [ALPHA Collaboration], JHEP 0107 
(2001) 048; 

M. Delia Morte, R. Frezzotti, J. Heitger and S. Sint [ALPHA Collabo- 
ration], JHEP 0110 (2001) 041. 

[2] Ph. Boucaud et al. [ETM Collaboration], Phys. Lett. B650 (2007) 304. 

[3] P. Dimopoulos et al, larXiv:0810.2873li [hep-lat]; P. Dimopoulos et al, 
PoS(LATTICE 2007)102; 
C. Urbach, PoS(LATTICE 2007)022; 
C. Michael, PoS (LATTICE 2007)122. 

[4] R. Frezzotti, G. Martinelli, M. Papinutto and G.C. Rossi, JHEP 0604 
(2006) 038 and PoS LAT2005:285, 2006. 

[5] B. Blossier et al. [ETM Collaboration], JHEP 0804 (2008) 020. 

[6] C. Alexandrou et al. [ETM Collaboration], Phys. Rev. D78 (2008) 
014509. 

[7] Ph. Boucaud et al. [ETM Collaboration], Comput. Phys. Commun. 179 
(2008) 695. 

[8] K. Cichy, J. Gonzalez Lopez, K. Jansen, A. Kujawa and A. Shindler, 
Nucl. Phys. B800 (2008) 94. 



33 



[9] R. Frezzotti, V. Lubicz and S. Simula [ETM Collaboration], 
larXiv:0812.4042l [hep-lat]. 

[10] B. Blossier et al. [ETM Collaboration], JHEP 0907 (2009) 043. 

[11] R. Frezzotti and G.C. Rossi, JHEP 0408 (2004) 007. 

[12] L. Scorzato, Eur. Phys. J. C37 (2004) 445. 

[13] S.R. Sharpe and J.M.S. Wu, Phys. Rev. D70 (2004) 094029 and ibidem 
D71 (2005) 074501. 

[14] G. Minister, JHEP 0409 (2004) 035. 

[15] F. Farchioni et al. [ETM Collaboration], Eur. Phys. J. C39 (2005) 421. 

[16] F. Farchioni, K. Jansen, I. Montvay, E.E. Scholz, L. Scorzato, A. 
Shindler, N. Ukita, C. Urbach and I. Wetzorke, Eur. Phys. J. C42 (2005) 
73. 

[17] F. Farchioni, K. Jansen, I. Montvay, E.E. Scholz, L. Scorzato, A. 
Shindler, N. Ukita, C. Urbach, U. Wenger and I. Wetzorke, Phys. Lett. 
B624 (2005) 324. 

[18] K. Jansen et al. [XLF Collaboration], Phys. Lett. B624 (2005) 334. 

[19] K. Symanzik, Nucl. Phys. B226 (1983) 187 and 205. 

[20] M. Liischer, S. Sint, R. Sommer and P. Weisz, Nucl. Phys. B478 (1996) 
365. 

[21] R. Frezzotti and G.C. Rossi, PoS(LATTICE 2007)277. 

[22] G. Minister and C. Schmidt, Europhys. Lett. 66 (2004) 652; 
G. Miinster, JHEP 0409 (2004) 035. 

[23] R. Frezzotti, talk given at the Ringberg Castle Workshop (April 2006), 
unpublished. 

[24] C.R. Allton et a/., Phys. Lett. B453 (1999) 30. 

[25] R. Babich, N. Garron, C. Hoelbling, J. Howard, L. Lellouch and C. 
Rebbi, Phys. Rev. D74 (2006) 073009. 

[26] C. Allton et al. [RBC-UKQCD Collaboration], Phys. Rev. D 78 (2008) 
114509; 

C. Aubin, J. Laiho and R.S. Van de Water. larXTv:0905.3947l [hep-lat]. 

34 



[27] P. Dimopoulos et al. [ETM Collaboration] PoS LATTICE2008:271,2008. 

[28] S. Coleman, J. Wess and B. Zumino, Phys. Rev. 177 (1969) 2239; 

C.G. Callan, S. Coleman, J. Wess and B. Zumino, Phys. Rev. 177 (1969) 
2247; 

S. Weinberg, Physica 96 A (1979) 327. 

[29] S. Weinberg, "The Quantum Theory of Fields" , Vol. 1, Cambridge Uni- 
versity Press (Cambridge, 1998). 

[30] J.F. Donoghue, E. Golowich and B.R. Holstein, "Dynamics of the Stan- 
dard Model\ Cambridge Univiversity Press (Cambridge, 1992). 

[31] M. Gell-Mann, R.J. Oakes and B. Renner, Phys. Rev. 175 (1968) 2195. 

[32] G. Martinelli, C. Pittori, C.T. Sachrajda, M. Testa and A. Vladikas, 
Nucl. Phys. B445 (1995) 81. 

[33] P. Dimopoulos, R. Frezzotti, G. Herdoiza, A. Vladikas, V. Lubicz, S. 
Simula and M. Papinutto, PoS LAT2007:241,2007. 

[34] C. Amsler et al. (Particle Data Group) PLB 667, 1 (2008). 

[35] B. Sheikholeslami and R. Wohlert, Nucl. Phys. B259 (1985) 572. 

[36] R. Frezzotti and G.C. Rossi, JHEP 0410 (2004) 070. 

[37] R. Frezzotti and G.C. Rossi, Nucl. Phys. B (Proc. Suppl.) 129 (2004) 
880. 

[38] S.R. Sharpe, Phys. Rev. D72 (2005) 074510. 

[39] M. Liischer, Commun. Math. Phys. 104 (1986) 177; Commun. Math. 
Phys. 105 (1986) 153 and Nucl. Phys. B354 (1991) 531. 

[40] M. Liischer, Nucl. Phys. B364 (1991) 237. 

[41] C.J.D. Lin, G. Martinelli, C.T. Sachrajda and M. Testa, Nucl. Phys. 
B619 (2001) 467. 

[42] S. Aoki and O. Bar, Phys. Rev. D70 (2004) 116011. 

[43] A.M. Abdel-Rehim, R. Lewis and R.M. Woloshyn, Phys. Rev. D71 
(2005) 094505. 



35 



[44] K. Jansen, M. Papinutto, A. Shindler, C. Urbach and I. Wetzorke [XLF 
Collaboration], Phys. Lett. B619 (2005) 184 and JHEP 0509 (2005) 
071. 

[45] J. Gasser and U.G. Meissner, Phys. Lett. B258 (1991) 219. 



36 



obs. O 


/? 


ajj q 


R s [0) 


m PS 


3.90 


0.0040 


0.185(44) 


m PS 


3.90 


0.0085 


0.139(51) 


m PS 


4.05 


0.0030 


0.120(110) 


m PS 


4.05 


0.0060 


0.120(42) 


/ps 


3.90 


0.0040 


0.04(6) 


/ps 


3.90 


0.0085 


-0.09(8) 


/ps 


4.05 


0.0030 


-0.03(6) 


/ps 


4.05 


0.0060 


0.01(5) 


my 


3.90 


0.0040 


0.022(68) 


my 


3.90 


0.0085 


0.021(44) 


my 


4.05 


0.0030 


-0.104(108) 


my 


4.05 


0.0060 


0.003(49) 


myfy 


3.90 


0.0040 


-0.07(18) 


myfy 


3.90 


0.0085 


-0.01(11) 


myfy 


4.05 


0.0030 


-0.31(29) 


Ulyfy 


4.05 


0.0060 


-0.12(13) 


m A 


3.90 


0.0040 


0.022(29) 


mA 


3.90 


0.0060 


0.001(21) 


m A 


3.90 


0.0085 


0.005(19) 


m A 


3.90 


0.0100 


0.007(19) 


m A 


4.05 


0.0030 


-0.004(45) 


m A 


4.05 


0.0060 


0.004(17) 


m A 


4.05 


0.0080 


0.000(18) 


m A 


4.05 


0.0120 


0.007(16) 



Table 1: The ratio -R S [C] of eq. (13. ip for the observables O indicated in the 
first column at the simulation parameters specified in the second and third 
column (lattice size is about L = 2.1 fm). The large statistical error on R s 
for mesonic observables at (3 = 4.05, a\i q = 0.0030 is due to limited statistics 
for the quark-disconnected contributions to the neutral meson correlators. 
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{P,a^ q )\L/a 


(3.9, .0040)| 24 


(3.9, .0064) | 24 


(3.9, .0085) | 2 4 


(3.9, .0040) | 32 


# meas. 


880 


240 


248 


184 


am^±\ L 


0.1359(7) 


0.1694(4) 


0.1940(5) 


0.1338(3) 


a 2 G n ±\ L V2 


0.1501(25) 


0.1581(16) 


0.1643(15) 


0.1484(11) 




0.103(4) 


0.134(10) 


0.163(8) 


0.107(7) 


am a \ L 


0.234(25) 


0.391(46) 


0.449(47) 


0.284(47) 


a 2 CiV2 


-0.021(4) 


-0.018(6) 


-0.011(6) 


-0.019(5) 




-0.12(2) 


-0.22(5) 


-0.25(6) 


-0.15(3) 




-3.0(6) 


-2.6(9) 


-1.6(9) 


-2.7(7) 




-0.10(2) 


-0.14(5) 


-0.10(4) 


-0.08(2) 


(P,afi q )\ L / a 


(4.05,. 0030) | 32 


(4.05, .0060)| 32 






# meas. 


192 


188 






am n ± \l 


0.1038(6) 


0.1432(6) 






a 2 G n ±\ L V2 


0.0898(18) 


0.0972(12) 






am^L 


0.090(6) 


0.125(6) 






am a \ L 


0.231(37) 


0.232(55) 






a 2 c t V2 


-0.014(3) 


-0.010(6) 






a 2 d x \/2 


-0.10(2) 


-0.09(3) 






Cir^a~ l \^2 


-4.0(9) 


-2.9(1.7) 








-0.13(2) 


-0.16(9) 







Table 2: Quantities entering eqs. ( 1C.46l) -( fC.49l) . The labels for the parameter 
sets have self-explanatory names and the corresponding columns follow the 
same order as the list of parameter sets in eqs. (1C50[) (1C.52[) . The statistical 
errors come from a standard bootstrap analysis on blocked data, so as to 
properly take into account autocorrelation of consecutive measurements. 
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